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ABSTRACT 

We investigate the local correlation betwen the 1.4 GHz radio continuum and 60 /im 
far-infrared (FIR) emission within the Large Magellanic Cloud (LMC) on spatial scales 
between 0.05 and 1.5 kpc. On scales below ~1 kpc, the radio-FIR correlation is clearly- 
better than the correlation of the cold gas tracers with either the radio or the FIR 
[ emission. For the LMC as a whole, there is a tight correlation between the radio and 

FIR emission on spatial scales above ~50 pc. By decomposing the radio emission into 
thermal and non-thermal components, however, we show that the scale on which the 
radio-FIR correlation breaks down is inversely proportional to the thermal fraction 
Q ■ of the radio emission: regions that show a strong correlation to very small scales are 

S—{ | the same regions where the thermal fraction of the radio emission is high. Contrary 

c/3 . to previous studies of the local radio-FIR correlation in the LMC, we show that the 

slope of the relation between the radio and FIR emission is non-linear. In bright star- 
forming regions, the radio emission increases faster than linearly with respect to the 
FIR emission (power-law slope of ~ 1.2), whereas a flatter slope of ~ 0.6 — 0.9 applies 
' more generally across the LMC. Our results are consistent with a scenario in which 

the UV photons and cosmic rays in the LMC have a common origin in massive star 
formation, but the cosmic rays are able to diffuse away from their production sites. 
Our results do not provide direct evidence for coupling between the magnetic field 
and the local gas density, but we note that synchrotron emission may not be a good 
tracer of the magnetic field if cosmic rays can readily escape the LMC. 

Key words: galaxies: individual (Large Magellanic Cloud) - ISM: galaxies - radio 
continuum: galaxies - infrared: galaxies - Magellanic Clouds 
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1 INTRODUCTION 

The tight, almost ubiquitous correlation between the 
far-infrared (FIR) and radio continuum emission in star- 
forming galaxies remains one of the most robust, yet 
puzzling, relatio nships in extragalactic astronomy (see e.g. 
lYun et ail 1200 ll see also Fig. 0. The correlation is essen- 
tially linear over five orders of magnitudes in luminosity 
with a dispersio n of less than 50%. I t applies to galaxies 
beyond z = 1 (lAppleton et all 12004) and encompasses a 
diverse range of galaxy types, including normal barred and 
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unbarred spirals, dwarf and irregular galaxies, st arbursts, 
Seyfe rts and radio-quiet quasars (for a review, see ICondonl 

The conventional explanation for the radio-FIR cor- 
relation invokes massive star formation. In this model, 
the FIR emission is mostly thermal emission from dust 
that has been heated by the ultra-violet (UV) radiation 
from young massive stars. These same stars power the 
thermal radio emission of star-forming regions, and rapidly 
evolve to the supernovae whose remnants are responsible 
for accelerating the cosmic-ray electrons that produce the 
non-thermal component of the radio emission. The radio 
emission of normal spiral galaxies at centimetre wavelengths 
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is predominantly non-thermal JCondonl Il99^) . while the 
total FIR emission comprises both a warm dust component, 
heated in Hn regions by massive (>20 M©) ionising stars, 
and a significant cool dust compo nent that is heated by 
the i nterstellar radiation field (e.g. IWalterbos fc Schwerind 
Il987f) . Efforts to decompose the radio emission into thermal 
and non-thermal components, and the FIR emission into 
warm and cool dust components, indicate that the overall 
radio- FIR correlation is constituted by a strong correlation 
between the thermal radio and the warm dust emission, 
and a weaker correlation between the non-the rmal radio 
and the cool dust emission (iHoernes et alJll998ft . 

Although a correlation between warm dust and thermal 
radio emission from ionised gas is hardly surprising, the 
mechanism by which the dust grains couple to synchrotron 
emission arising through the interactions of relativistic 
electrons with the ambient ma gnetic field remains highly 
uncertain. iHelou fc Bicavl lll993t) have suggested a coupling 
between magnetic field strength and gas density, which 
would ensure a tight radio-FIR correlation in spite of the 
sensitive dependence of the sychrotron emission on the 
m agnetic field. Variatio ns o n this idea have bee n discussed 
bv lHoernes et alJ il99ST) and lMurgia et al.l J2005h . the latter 
also connecting the radio-FIR correlation to the observed 
CO-radio correlation, as CO traces the dense molecular gas. 
An important constraint on these models is the scale on 
which the correlation breaks down, as this would shed light 
on the physical mechanism of the coupling (e.g., static pres- 
sure equilibrium predicts that the correlation should hold 
on scales greater than the scale height of a galaxy's gas disc). 

Observationally, FIR and radio maps of nearby spiral 
galaxies generally show good correspondence, although 
studies of nearby spir als are limited to scales above a few 
hund red parsecs (e.g. IHoernes et al.lll993 iHippelein et alJ 
2003). On much smaller scales, detailed maps of the gas 
within several kiloparsecs of the Sun indicate that the 
radio emission from star-forming regions is predominantly 
thermal, and does not follow the radio-FI R correlation for 
integ r ated galaxy luminosities (e.g. iBoulanger fc Peraultl 
ll98Sl: lHaslam fc Osborne! fl987^l . Galactic supernova rem- 
nants do not seem to be well-correlated with FIR or 
mid-infr ared (MIR) emission, or with sites of star for - 
mation dCohen fc Greenl l200ll : IWhiteoak fc Greer] Il99rl) . 
A tight, linear correlation between the integrated radio 
and FIR emission in galaxies is thus not at all expected 
from observations of the local interstellar medium (ISM). 
In light of this conflict, and the different effective scales 
of proposed mechanisms for the origin of the radio-FIR 
correlation, it seems timely to study the radio and FIR 
emission of a galaxy on intermediate spatial scales between 
those accessible in the Milky Way and in nearby spirals. 

The aim of this study is to investigate the radio-FIR 
correlation within the Large Magellanic Cloud (LMC). 
The LMC is a gas-rich dwarf galaxy that exhibits clear 
signs of active star formation. Our global view of the 
LMC is complemented by sufficient spatial resolution to 
study the ISM structure in detail over a wide ra nge of 
wave bands. At an assumed distance of 50.1 kpc (lAlvesI 
12004) . 1 arcmin corresponds to 14.5 pc. The proximity of 



the LMC thus allows a detailed comparison of the radio and 
FIR emission on scales down to 20 pc (for the highest res- 
olution radio and FIR datasets that are currently available) . 

In addition to standard pixel-by-pixel techniques, we 
use a wavelet cross-correlation method for our analysis 
that allows us to quantify the correlation between different 
gas and dust tracers on spatial scales between ~0.04 
and 2.0 kpc. To date, wavelets have enjoyed a limited 
(though rapidly increasing) application within astrophysics, 
but their ability to isolate signal simultaneously in the 
Fourier and spatial domains makes them ideal for analysing 
scale-dependent correlations in astronomical images. 

The structure of this paper is as follows. In Section 
we describe our data sets and reduction procedures. Our 
analysis methods and results are presented in Section [3] We 
discuss our results in the context of current models for the 
origin of the radio-FIR correlation in Section^] Our conclu- 
sions are presented in Sectional 



2 OBSERVATIONAL DATA 

The images we have used for this study are described in 
detail below, and summarized in Table Q Originating from 
a variety of instruments, they cover different fields of view 
at different angular resolutions and pixel sizes. In order to 
compare the data, all maps were regridded to a common 
pixel size of 20 arcsec and a common J2000 map centre 
of a=5 h 20 m 38 , !6, <5=-68°41'22". Images of 1350 2 pixels, 
covering a roughly 7? 5 x 7? 5 region, were produced using 
an NCP projection, which is essentially an orthographic 
(SIN) projection with a tangent point at the north celestial 
pole. The final regridded, truncated images are shown 
in Fig. [5] Note that the region thus defined, which we 
refer to as the "whole LMC," is completely covered by 
all maps except for the CO map, which has a somewhat 
irregular sky coverage (Fig. |5J . Total fluxes and flux densi- 
ties measured within the selected region are given in Table0 

Since we are interested in spatial variation of the 
radio-FIR correlation across the LMC, the maps were 
divided into 16 sub-regions, as shown in Fig. Each 
sub-region is 1?35 x 1?35 in size, with 256 x 256 pixels and 
a pixel scale of 19 arcsec. The centre positions of the 16 
sub-region maps are listed in Table |2] along with brief notes 
about the individual fields. Note the wide variation in flux 
density among sub-regions in the 1.4 GHz and 60 /xm maps: 
for instance, sub-region 8, which covers the 30 Doradus 
region, contains roughly a third of the total radio and FIR 
emission in the LMC. 



2.1 Far-Infrared Data 

Infrared Astronomical Satellite (IRAS) mosaics of the LMC 
in the 60 /an and 100 /jm bands were kindly provided to 
us by M. Braun and have been discussed by iBraun et alJ 
(1998). These data had been processed using the HIRAS 
algorithm, which uses the Pyramid Maximum Entropy 
(PME) method to improve upon the nominal resolution 
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Figu re 1. Plot o f 1.4 GHz vs IRAS 60 /im flux density (left), and 1.4 GHz radio luminosity vs. IRAS 60 (im luminosity (right) for the 
lYun et alj 1200 ll) sample of 1809 galaxies. In both panels, the dashed line is a linear least squares fit to the data. The LMC is indicated 
as an open square. 



Table 1. Properties of the datasets used in this study, prior to 
any convolution or regridding. The fluxes and flux densities are 
measured over a roughly 7? 5 X 7? 5 region covering the main body 
of the LMC. 



Band 


Telescope(s) 


Resolution 


Pixel scale 


Flux Density 






(arcmin) 


(arcmin) 


(Jy) 


60 fim 


IRAS (HIRAS) 


1.3 


0.25 


9.64 x 10 4 


60 fim 


IRAS (IRIS) 


4.Q 


1.5 


9.61 x 10 4 


100 /im 


IRAS (HIRAS) 


3.0 


0.25 


2.09 x 10 s 


100 /^m 


IRAS (IRIS) 


4.3 


1.5 


2.04 x 10 5 


6.3 cm 


Parkes 


4.8 


1.8 


2.96 x 10 2 


21 cm 


Parkes 


15.2 


5.0 


4.71 x 10 2 


21cm 


ATCA+Parkes 


0.67 


0.22 


4.26 x 10 2 


Band 


Telescope(s) 


Resolution 


Pixel scale 


Flux 






(arcmin) 


(arcmin) 


(Jy km s -1 ) 


Hi 


ATCA+Parkes 


1.0 


0.33 


3.98 X 10 s 


CO 


NANTEN 


2.6 


2.0 


8.80 x 10 5 



can be found on the IRAS Documentation website. 1 

To verify the flux density calibration of the HIRAS 
data, IRIS 60 fim and 100 fim maps were obtained from Dr. 
Mivil le-Deschenes at CITA iMiville-Deschenes fc Laeachd 
2005). The IRIS images are at the same resolution as the 
original IRAS survey data, but have better zodiacal light 
subtraction and destriping than the equivalent IRAS Sky 
Survey images. The flux density calibration of IRIS is con- 
sistent with DIRBE data, but we found it was necessary to 
apply a small zero-level correction to the maps. Measuring 
the average brightness for regions of sky away from the 
LMC, we determined a background level of 0.85 MJy sr _1 
for the 60 /im map and 3.0 MJy sr _1 for the 100 /an 
map. We subtracted these constant values from the maps. 
It is worth noting that the background levels measured 
for the IRIS maps are very similar to the background 
levels measured for the HIRAS maps, which increases our 
confidence in the flux density calibration of the HIRAS 
maps. A complete description of IRIS data products can be 
found on the IRIS website. 2 



of the original IRAS survey data (4.0T0.2 arcmin for 
the 60 /mi band, an d 4.3T0.2 for the 100 fim band) 
l|Bontekoe et alJll994lL The PME algorithm relies on the 
local sampling pattern in the IRAS survey, so the final 
resolution of HIRAS images varies spatially across a map. 
By fitting Gaussian profiles to point-like sources, we found 
that the resolution of the HIRAS maps was ~1.3 arcmin 
for the 60 fim data and ~3.0 arcmin for the 100 fim data. 
Measuring the average brightness for regions of sky away 
from the LMC, we determined a background level of 1.25 
MJy sr _1 for the 60 fim map and 3.0 MJy sr _1 for the 
100 fim map. These constant values were subtracted from 
the maps. A complete description of IRAS data products 



2.2 Radio Continuum Data 

We used a map of the 1.4 GHz radio continuum emission 
from the LMC produced by combining data from the ATNF 
Parkes 64 m telescope and the 6-element Australia Tele- 
scope Compact Array (ATCA). 3 Details of the o bservations 
and data reduction were presented in papers bv lKim et alJ 

1 http:/ /irsa. ipac.caltech.edu/IRASdocs/toc. html 

2 http:/ /www. cita.utoronto.ca/~mamd/IRIS/IrisOverview. html 

3 The Australia Telescope is funded by the Commonwealth of 
Australia for operation as a National Facility managed by CSIRO. 
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Figure 2. Maps of the Large Magellanic Cloud at different wavelengths, overlaid by the sub-region fields described in the text. The 
images correspond to HIRAS 60 (im emission (top left), HIRAS 100 fim emission (top right), ATCA+Parkcs 1.4 GHz radio continuum 
(middle left), median filtered ATCA+Parkes 1.4 GHz radio continuum (middle right), ATCA+Parkes Hi integrated intensity (bottom 
left), and NANTEN CO J=l-0 integrated intensity (bottom right). A square-root intensity scale has been used in displaying all images. 
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Table 2. Position centres, absolute and relative fluxes, and absolute and relative flux densities of the 16 sub-regions. Flux densities 
for the Parkes 1.4 and 4.8 GHz images are presented in Jy. The IRIS 60 (im and 100 fim flux densities are presented in 10 3 Jy. The 
fluxes for the Hi and CO are presented in 10 4 Jy km s — 1 . The relative contribution for a sub-region (shown in parentheses) is given as 
a percentage of the total emission in the entire region shown in Fig. [5] 



Reg. 


R. A. 


Dec. 


Si, 4 








/hi 


^CO 


Notes 




(J2000) 


(J2000) 


Jy 


Jy 


10 3 Jy 


10 3 Jy 


10 4 Jy km s- 1 


10 4 Jy km s- 1 




1 


4:52:56 


-70:49:30 


4.4 (1%) 


1.4 (0%) 


0.6 (1%) 


2.7 (1%) 


0.9 (2%) 


0.6 (1%) 




2 


5:09:22 


-70:49:30 


9.0 (2%) 


4.1 (1%) 


1.4 (1%) 


4.2 (2%) 


0.9 (2%) 


3.4 (4%) 




3 


5:25:49 


-70:49:30 


20.7 (4%) 


12.3 (4%) 


2.7 (3%) 


6.8 (3%) 


1.5 (4%) 


6.6 (7%) 




1 


5:42:15 


-70:49:30 


25.4 (5%) 


13.7 (5%) 


2.9 (3%) 


10.1 (5%) 


4.1 (10%) 


12.1 (14%) 




5 


4:54:29 


-69:28:30 


14.0 (3%) 


9.2 (3%) 


3.8 (4%) 


9.1 (4%) 


1.5 (4%) 


5.0 (6%) 




6 


5:09:53 


-69:28:30 


21.8 (5%) 


13.3 (5%) 


6.0 (6%) 


12.2 (6%) 


1.5 (4%) 


8.1 (9%) 




7 


5:25:18 


-69:28:30 


39.4 (8%) 


27.0 (9%) 


8.1 (8%) 


16.0 (8%) 


1.6 (4%) 


6.3 (7%) 


N132D 


8 


5:40:42 


-69:28:30 


129.6 (28%) 


107.2 (36%) 


36.9 (38%) 


53.5 (26%) 


4.9 (12%) 


16.5 (19%) 


30 Dor 


9 


4:55:51 


-68:07:30 


8.3 (2%) 


4.7 (2%) 


2.0 (2%) 


5.4 (3%) 


1.6 (4%) 


3.9 (4%) 




10 


5:10:21 


-68:07:30 


21.0 (4%) 


11.5 (4%) 


3.5 (4%) 


8.5 (4%) 


2.1 (5%) 


4.7 (5%) 




11 


5:24:50 


-68:07:30 


34.2 (7%) 


24.0 (8%) 


7.6 (8%) 


15.9 (8%) 


2.6 (7%) 


8.9 (10%) 




12 


5:39:20 


-68:07:30 


24.7 (5%) 


19.0 (6%) 


4.8 (5%) 


11.2 (5%) 


2.3 (6%) 


3.5 (4%) 




13 


4:57:03 


-66:46:29 


14.1 (3%) 


9.8 (3%) 


3.7 (4%) 


8.2 (4%) 


1.9 (5%) 


5.3 (6%) 


Nil 


14 


5:10:45 


-66:46:29 


11.8 (3%) 


4.8 (2%) 


1.6 (2%) 


4.6 (2%) 


1.7 (4%) 


2.3 (3%) 


PKS 0515-674 


15 


5:24:26 


-66:46:29 


17.8 (4%) 


10.6 (4%) 


2.6 (3%) 


6.6 (3%) 


1.5 (4%) 


2.9 (3%) 




16 


5:38:08 


-66:46:29 


13.1 (3%) 


8.7 (3%) 


1.9 (2%) 


5.0 (2%) 


1.0 (3%) 


0.7 (1%) 





lll998t I2003T) . The angular resolution of the image is 40 
arcsec, corresponding to a spatial resolution of 10 pc. The 
image suffers from ring-like artefacts at the ~0.5% level, 
presumably due to errors in the primary beam model used 
during the mosaicking and deconvolution process. These er- 
rors become significant close to bright compact sources such 
as the 30 Doradus region on the eastern side of the galaxy, 
limiting the sensitivity that can be achieved in these regions. 

The 1.4 GHz image contains a large num ber of point 
sourc es, of which >90% are background AGN iMarx et alJ 
Il997l . Filipovic et al., in prep.). For the purpose of our 
study, it was desirable to remove these, since they lead to 
a spurious decorrelation on small scales. We found that 
most of the background radio sources could be eliminated 
by median filtering the image. The filtering operation, 
implemented in the GIPSY routine mfilter, moves a 
2^5x2.' 5 window across a map, replacing the central pixel 
value (S C pix) with the median value of the window (Smed) if 
\S C pix — S m ed\ > S me d + 1 mjy The 1 mjy offset prevents 
unnecessary filtering in noisy regions where the median is 
close to zero. The final median filtered radio continuum 
map is shown in Fig. [5] (middle right). For the whole 
LMC, the filter removed ~10% of the total emission in the 
1.4 GHz map. For individual sub-regions, the fraction of the 
total emission removed by median filtering varies between 
a few percent, for sub-regions with strong emission from 
the LMC itself, to ~40% in sub-region 1. The ring artefacts 
surrounding extremely bright sources such as 30 Doradus 
(subregion 8), N132D (subregion 7) and PKS0515-674 
(subregion 14) were not removed by median filtering. 

Lower resolution radio maps of the LMC at 1.4 
and 4.8 GHz were obtain ed from the Parkes surveys of 
lHavnes et alJ lll986t ll99lT) in order to estimate the large 
scale variation in the thermal fraction of the radio emission. 



The angular resolutions of the single-dish 1.4 and 4.8 GHz 
radio maps are 15.2 and 4.8 arcmi n respectively. Furthe r 
details are presented in the paper bv lFilipovic et al] lll995l) . 

2.3 Hi Data 

iKim et all <l2003l) have combined Parkes single-dish neutral 
hydrogen observations of the LMC with an ATCA aper- 
ture synthesis mosaic to generate an Hi data cube for the 
LMC with limiting angular resolution of 1 arcmin. The final 
dataset is sensitive to spatial scales above 15 pc across the 
heliocentric velocity range 190 km s" 1 to 388 km s _1 . The 
velocity resolution of the Hi data is 1.65 km s~ . An inte- 
grated intensity map was produced from the combined data 
cube by integrating over the heliocentric velocity range of 
195 to 355 km s _1 and clipping values with absolute value 
< 2cr. The use of clipping does bias the resulting map against 
low-level emission below the clipping value, but in total the 
excluded flux is only a small fraction (~8%) of the flux in 
an undipped integrated map. 

2.4 Molecular Gas Data 

A survey in 12 CO(J=1-0) has been carried out by NAN- 
TEN, a 4m radio telescope operated by Nagoya University 
at Las Campanas Observatory i n Chi le. This new sur- 
vey, introduced by iFukui et alJ i200ll) . has achieved a 
sensitivity equivalent to N(H.2) ~ 1.0 x 10 21 cm" 2 , a 
facto r of ^3 better than the original survey llFukui et alJ 
Il999t iMizuno et all 1200 ill . The data covers an irregularly 
shaped area of ~ 6° x 6° centered on the B1950 position 
(5 h 20 m 00 s , -69°00'00"), and has revealed a distribution 
of 259 giant molecular clouds with 27,000 observed posi- 
tions iFukui et al.ll200ll) . The half-power beam width of 
the telescope was 2.6 arcmin and data were taken with 
a 2 arcmin grid spacing. The spectrometers were two 
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acousto-optical spectrometers with 2048 channels each. The 
narrow-band (NB) spectrometer has a velocity coverage 
of 100 km s _1 and a resolution of 0.1 km s _1 , whereas 
the wide-band (WB) spectrometer has a velocity coverage 
of 650 km s _ and a resolution of 0.65 km s _1 . Out of 
the 27,000 positions, 6229 were observed with the NB 
spectrometer, while the rest were observed with the newly 
developed WB spectrometer. 

An integrated intensity map was produced from the 
NANTEN data cube as follows. The velocity channels were 
first binned to 1.95 km s^ 1 . To isolate regions of emission, 
a blanking mask was defined using the 3-<r contour of a 
smoothed version of the data cube; this cube had been 
generated by convolution of the binned cube with a Gaus- 
sian kernel with FWHM=5.'2. To create the integrated 
intensity map, the mask was applied to the original binned 
cube, which was then summed over local-standard-of-rest 
velocities of 200 to 300 km s" 1 . 



3 ANALYSIS 

For this study, we are interested in both the structural cor- 
respondence between images taken in different wavebands 
and changes in the ratio of radio to FIR emission across the 
LMC. These considerations led us to adopt two techniques 
for our cross-correlation analysis: a wavelet method, which 
allows us to probe morphological correlations across the 
whole range of spatial scales present in each image, and 
a more traditional pixel-by-pixel comparison, which is 
sensitive to changes in the emission ratio as a function of 
position. We examine the FIR/radio ratio across the LMC 
in order to compare our results with other studies of the 
radio-FIR correlation in nearby disc galaxies. 

Prior to conducting our analysis, maps of the whole 
LMC and the individual sub-regions were smoothed to 
a common resolution. For comparisons between the H I 
and 1.4 GHz continuum data, the continuum data were 
smoothed to 1 arcmin, i.e. the resolution of the Hi data. 
For comparisons with the CO data, maps were smoothed to 
3 arcmin, slightly larger than the FWHM of the NANTEN 
beam since the map was not Nyquist sampled. Comparisons 
involving the 60 fim and 100 /im maps were less straight- 
forward, due to the spatially varying HIRAS resolution in 
response to survey coverage and signal-to-noise. We tested 
for potential resolution-dependent effects on our results 
by repeating our cross-correlation analysis with a HIRAS 
map smoothed to a range of resolutions between 1.5 and 5 
arcmin. This simple validation test indicated that for most 
sub-regions the cross-correlation spectrum at small scales 
can be sensitive to changes in resolution, but only when the 
variation in the resolution is quite large (i.e. a factor of ~2). 
Visual inspection of the HIRAS data, however, suggested 
that the resolution of the map around 30 Doradus may be 
significantly lower than our nominal value of 1.3 arcmin. 
Further discussion of possible resolution effects on the 
cross-correlation results for this sub-region is presented in 
Section l3~4l 



All the data preparation and analysis tasks in this study 
were performed with the Miriad, Gipsy and Yorick soft- 
ware packag es llSault et alJll995t Ivan der Hulst et al.lll992l 
lMunrolll9951) . 4 



3.1 Pixel- by-pixel Analysis 

The simplest measure of the correlation between two images, 
fi(x,y) and fi{x,y), with the same angular resolution and 
the same number of pixels is Pearson's linear correlation 



at each pixel: 



£(A* - (A))(/2* - {h)) 
- </i» 2 - </ 2 » 2 



(1) 



Identical images should have r p = 1; images that are per- 
fectly anti-correlated should have r v = — 1. The formal error 
on the correlation coefficient depends on the strength of the 
correlation and the number of independent pixels, n, in an 



image: 



(2) 



In our case, n = 2601 for the sub-region maps and 
n — 72900 for the whole LMC, so the formal error is small. 
The true error in the correlation coefficient is dominated 
by systematic errors in the data, such as calibration and 
zero-level uncertainties, which are not so easily quantified. 

The correlation coefficients for the median filtered 
1.4 GHz and 60 /im emission in the whole LMC and 
individual sub-regions are presented in Table |H] The cor- 
relation between log Si. 4 and log S60 was calculated, with 
flux densities expressed in Janskys. To ensure statistical 
independence, one in every five pixels along both axes were 
used for the calculation of r v , corresponding to a separation 
of 1.7 arcmin for the whole LMC map, and 1.6 arcmin for 
the sub-region maps. For the whole LMC, r p =0.86, while 
the median of the correlation coefficients for the sub-regions 
is 0.73. The correlation coefficients of individual sub-regions 
show a large scatter about this value, ranging between 0.36 
and 0.91. As measured by Pearson's correlation coefficient, 
however, the correlation between the 1.4 GHz and 60 /im 
(or 100 fira) emission across the LMC is much stronger than 
the correlation found for other image pairs (see Table 0). 

Given that a strong correlation is present, we can 
consider how the slope of the correlation varies between the 
individual sub-regions and the LMC as a whole. Scatter 
plots in logarithmic coordinates showing the flux density 
of the 60 fim against the median filtered 1.4 GHz radio 
emission for the whole LMC and the individual sub-regions 
are shown in Figs. |3] and 0| respectively. As for Pearson's 
correlation coefficient, the plots are generated from maps 
that were smoothed to 1.5 arcmin resolution; the pixels that 
are used to construct the scatterplots are separated along 
both axes by 1.7 arcmin for the whole LMC maps, and 1.6 
arcmin for the sub-regions. As many data points can occupy 
the same region in a plot, we indicate the distribution of 

4 http://yorick.sourceforge.net 
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Figure 3. Pixel- by- pixel correlation between the 60 fim and me- 
dian filtered 1.4 GHz maps of the LMC. The solid line is the WLS 
fit to the data (slope = 1.18), while the dashed line is the OLS fit 
(slope = 0.75). The white contours indicate the density of data 
points, with contour levels 5, 30, 60 and 90% of the distribution's 
peak density. The black contours, also with contour levels at 5, 30, 
60 and 90% of the noise's peak density, indicate the distribution 
that would be expected for Gaussian noise with RMS appropriate 
to the 1.4 GHz and 60 (im images. The open square represents the 
average surface brightness for the whole LMC, calculated using 
the total ATCA+Parkes 1.4 GHz and IRIS 60 fim flux densities 
listed in Table HI and a total galaxy area of 7? 5 X 7? 5. 

the data points by binning the data into a 2-D histogram 
and showing contours of the 2-D distribution. The mesh 
used to define the contour levels consists of 60 x 60 cells, 
evenly distributed in logarithmic space between 10 -7 and 
10 3 Jy arcmin -2 for the FIR emission, and 10~ 9 and 10 1 
Jy arcmin -2 for the radio emission. 

For both the individual sub-regions and the whole 
LMC, we find that the local correlation between the radio 
and the 60 /im emission is reasonably good. We used two 
methods for calculating a linear fit to the pixel-by-pixel 
plots: weighted least squares (WLS) and the ordinary 
least squares (OLS) bisector. The WLS fit (solid line), 
implemented usin g the Numerical Recipes fitexy routine 
(IPress et al.lll992T) . assumes an uncertainty of 2.3 x 10 -4 Jy 
arcmin -2 for each 1.4 GHz measurement, and 2.7 x 10 -2 Jy 
arcmin -2 for each 60 /im measurement. These uncertainities 
were estimated from the RMS of a blank region of sky in 
the 1.4 GHz and 60 /im LMC maps at 1.5 arcmin resolution. 
In logarithmic space, the high intensity pixels are thus 
given greater weighting due to their smaller relative errors 
(note that this assumes the errors are mainly additive 
rather than multiplicative). The OLS bisector (dashed 
line), implemented usin g the slopes routine provided by 
iFeieelson fc Babul <ll992T) . assigns equal weight to each pixel 
and treats both variables symmetrically by determining the 
line which bisects the standard OLS solutions of Y on X and 



X on Y. Measurement uncertainties are assumed to be unim- 
portant compared to the intrinsic scatter in the data. We 
therefore expect that the WLS method will provide a better 
fit to the high-intensity pixels, whereas the OLS bisector 
will provide a better fit to the overall distribution of pixels. 
In order to prevent noise pixels contributing to our results, a 
3-cr cut was applied to the data prior to determining the fits. 

For the whole LMC, the WLS method determined a 
slope of 1.18 to the data, while the OLS bisector indicated 
a flatter slope (0.75). Similar values are found for sub- 
regions 8 and 13, which contain the 30 Doradus and Nil 
star-forming complexes respectively. This is not surprising, 
since sub-region 8 is responsible for 38% (28%) of the 
LMC's total 60 /J,m (1.4 GHz) emission, and contributes 
significantly to the slope of WLS fit that we determine for 
the galaxy as a whole. The slopes of the WLS and OLS 
bisector fits in other sub-regions (listed in Table |3] and 
plotted in Fig. show a range of values between 0.55 and 
1.29. The median slope of the WLS fits is 0.78, while the 
median slope of the OLS bisectors is 0.75. Contrary to a 
previous study of the radio-FIR correlation in the LMC 
ilXu et alj|l992T) . our results suggest that locally there are 
two independent correlations between the radio and FIR 
emission, both of which are non-linear. The correlation 
in the bright star-forming regions is steeper than a linear 
correlation (slope~l.l-1.3), and it is superimposed on 
a more general correlation that is flatter than a linear 
cor relation ( s lope^ 0.6-0.9). We compare our results to 
the IXu et all 1^992) result in Section 14.11 and discuss the 
variation of the WLS fit with the thermal fraction of the 
radio emission in Section f3. 51 

Besides sub-regions 8 and 13, which contain 30 Dor 
and Nil and hence have a significant population of high- 
intensity pixels, the WLS fit and OLS bisector show greatest 
divergence in sub-region 1, where a large fraction of the 
emission is close to the sensitivity limit and the fit is there- 
fore poorly constrained. The small congregation of points 
that show a high radio to FIR ratio in sub -region 1 is due to 
the b right supernova remnant 0450-70.9 (iMathewson et all 
Il985h . 

lYun et alJ (1200 ll) derived a relation between the inte- 
grated 1.4 GHz and 60 fim emission from star- forming galax- 
ies given by 

log(Li. 4 /W Hz -1 ) = (0.99±0.1)log(L6o/i©)+(12.07±0.08) .(3) 

We used the data in their Table 1 to re-derive this relation 
in flux density units; a linear least squares fit gives 

log(Si. 4 /mJy) = 0.98 log(S 60 /Jy) + 0.86 . (4) 

These relations are shown as dashed lines in Fig. ^ O n both 
panels , the LMC lies slightly below fits to the fYun et alJ 
(2001) data. Furthermore, except for sub-regions 8 and 13, 
the slopes that we fit to the data points are flatter than 
a linear correlation. This causes the FIR/radio ratio to 
decrease for lower brightness regio ns, as has bee n note d 
in four nearb y spira l galaxies by iMurphv et alJ feOOEl) . 
IMurphv et alJ ll2005l) attribute this non-linearity of the 
local radio-FIR correlation to the diffusion of cosmic-ray 
electrons (CRe - s) away from star- forming sites, i.e. an 
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Figure 4. The pixel-by-pixel correlation between the 60 fim and median filtered 1.4 GHz radio maps for individual sub-regions. The 
linear fits and noise contours (black) are the same as in Fig. |!|] The white contours represent the density of the data points. The contour 
levels are 5, 30, 60 and 90% of the peak density in sub-region 11. 



individual star-forming complex will show an infrared 
excess due to local dust heating by UV photons, but it 
also generates an extended non-thermal radio halo, whose 
radius depends on the mean age of the CRe - population. 



3.2 FIR/radio ratio map 

In order to determine whether variations in the radio-FIR 
correlation correspond to particular physical structures in 
the LMC, we constructed a logarithmic FIR/radio ratio map 
("g-map") from the IRIS 60 and 100 /im maps, and the 



ATCA+Parkes 1.4 GHz radio map. These input maps were 
each smoothed to a resolution of 4.3 arcmin, i.e. the natural 
resolution of the 1 00 fim image. Following the convention of 
iHelou et alJ il985t) . q is defined as 

q = l0g ( 3 .75x ™Wm-0 - l0g ( wm-^Hz-0 ' (5) 

For each pixel, the total FIR surface brightness was deter- 
mined from a combination of the IRIS 60 an d 100 fj,m maps 
according to Equation 6 of lYun et alJ (l200lT) . 

FIR= 1.26 x 10" 14 (2.58B 6 o + Bioo)Wm" 2 , (6) 
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Figure 5. Images of the logarithmic FIR/radio ratio, q, constructed from the combined IRIS 60 and 100 /an maps and a smoothed 
version of the ATCA+Parkes 1.4 GHz radio map. The left panel shows regions where q is less than its median value of 2.45, whereas the 
right panel shows regions where q > 2.45. Major structural features of the LMC are indicated. 



where Bqo and -Bioo are the 60 and 100 /im surface bright- 
ness in Jy pix -1 respectively. 

The resultant map is shown in Fig. [5] For clarity, the 
left panel displays pixels for which q > 2.45, and the right 
panel displays pixels where q < 2.45. Only pixels detected 
above the 3cr level in the radio map were included in the 
calculation, which produces an irregular-shaped mask in the 
final map. The median q across the whole LMC is 2.45±0.22, 
where 0.22 is the semi-interquartile range. For the LMC as 
a whole, q = 2.55, calculated using the total flux density of 
the IRIS 60/100^m and c ombined ATCA+ Parkes 1.4 GHz 
maps. For comparison, the lYun et alJ feOOlT) galaxy sample 
has a mean q — 2.34. Within the map boundaries where we 
calculate q, it thus appears that either the radio emission in 
the LMC is slightly unde rluminous c ompar ed to the value 
that we expect from the lYun et ail i200lT) correlation, or 
that the FIR emission is slightly overluminous (see Fig- El ■ 

In contrast to so me high-resolution g -maps of nearby 
disc galaxies (see e.g. iMurphv et al1l2005f) . the most strik- 
ing structures in the g-map of the LMC do not correspond 
to the brightest structures in the input FIR and radio maps 
(see Fig.|^J. However, inspection of the g-map shows that the 
variation in the FIR/radio ratio across the LMC is spatially 
organized to some extent. In particular, the stellar bar region 
and the western edge of the LMC show elevated FIR/radio 
ratios, while the eastern side of the LMC shows numerous 
patches of low q values. The supershell LM C4 in the north - 
east is the largest of these low q regions <MeaburrJll98fl . 
Both 30 Doradus and Nil have a pronounced ring-like mor- 
phology in the g-map. The local peak in the original radio 
and FIR maps is surrounded by a ri ng of enhanc e d q va lues. 
An earlier low-resolution study by iKlein et al. ( 1989J) has 
reported that LMC Hn regions have an average q = 2.78. 
Our data suggests that at higher angular resolution, the q 
value of an HII region can be separated into two distinct 



components: a compact central region where the FIR and 
(mostly thermal) radio emission are high but q is quite close 
to the mean q across the whole LMC, and a ring of enhanced 
q values that arises through the different scale-lengths of 
the FIR and radio emission associated with the Hn region. 
This morphology is consistent with the statistical detection 
of a "dip - and-r ing" structure around LMC Hn regions by 
IXu et all Jl992t . 



3.3 Wavelet Analysis 

Wavelet analysis involves the convolution of an image with 
a family of self-similar basis functions that depend on scale 
and location. It can be useful to consider the wavelet trans- 
form as the general case of the Fourier transform, where 
the oscillatory functions are localised in both time and fre- 
quency. The family of basis functions is generated by di- 
lating and translating a mother function, which is called 
the analysing wavelet. In this paper we make use of the 
two-dimensional continuous wavelet transform, which can 
be written as 

/oo p OO / i \ 

/ /(*>* dx' (7) 

- oo J — oo \ / 

Here /(x) is a two-dimensional function, such as an image, 
%[)* is the complex conjugate of the dilated and translated 
analysing wavelet, and w(a) is a normalising function. 
a is the scale size of the wavelet. Note that the term 
"continuous wavelet transform" is used in the literature 
to distinguish this transform from a "discrete" trans- 
form which employs an orthogonal set of basis functions; 
however, there is no fundamental difficulty applying the for- 
mer to discretely sampled data such as astronomical images. 

Our use of wa velet transforms follows closely the work of 
iFrick et all (l200lT) , who pioneered their application to galaxy 
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Figure 6. The wavelet (filled squares) and Fourier (solid line) 
power spectra for a 10" (FWHM) Gaussian. 



images. Our analysing wavelet is their "Pet Hat" which can 
be written in Fourier space as 



C(a) cos 2 (f log 2 2|k|a) 







i 

4 a 

|k| 



Ikl > 



(8) 



This function picks out an annulus in the Fourier plane 
centred at k=(2a)~ 1 , and provides both simplicity and good 
separation of scales. The normalising factor C(a) is chosen 
to provide unit total flux density in the wavelet image (or 
unit total flux, in the case of the Hland CO data). C(a) is 
evaluated numerically to be 1.065a/<5, where a and 8 are the 
scale size (defined below) and pixel size expressed in radians. 

Throughout this paper we adopt the convention that 
an angular scale size a in radians corresponds to a spatial 
frequency of (2a)" 1 wavelengths. This is based on the 
observation that a positive, Gaussian-like structure will 
tend to couple with half of a sine wave rather than a full 
period, however it differs from the usual convention by a 
factor of 2. As an illustrative example, we show in Fig.|S|the 
wavelet and Fourier spectra for a 10" (FWHM) Gaussian; 
the spatial frequency is shown on the top axis and the 
corresponding spatial scales are shown on the bottom axis. 
In Table [3] and the following sections, we express a in 
arcseconds, which are a more appropriate unit for the size 
scales of interest here. The Fourier spectrum is defined in 
the usual way as Ef(k) = |/(k)| 2 . 

Expressing the convolution of Equation Q in terms of 
Fourier transforms, the wavelet filtered image is calculated 



W(a,x) 



1 

AT2 



/(k)V*(ka) 

(ki,k 2 ) 



-2jri(k-x)/iV 



(9) 



Table 3. Summary of the wavelet and pixel-by-pixel correlations 
between the median filtered 1.4 GHz and the 60 £tm images for 
the whole LMC and individual sub-regions. Columns 2 and 3 list 
the angular (arcsec) and linear (parsec) spatial scales at which 
the wavelet cross-correlation falls below 0.75. Column 3 lists the 
Pearson's correlation coefficient, r p . Columns 4 and 5 list the 
WLS and OLS-bisector slopes. Column 6 lists the thermal fraction 
fth of the radio emission at 1.4 GHz, assuming a non-thermal 
spectral index of —0.7. 



Region log(a/arcsec) Scale (pc) 



slopey 



slopeoLS fth 



LMC 


2.3 


50 


0.86 


1.18 


0.75 


0.45 


1 


3.0 


>240 


0.36 


0.89 


1.17 


<0 


2 


2.7 


120 


0.57 


0.60 


0.64 


0.07 


3 


2.4 


60 


0.79 


0.76 


0.70 


0.37 


1 


3.0 


>240 


0.62 


0.64 


0.68 


0.26 


5 


2.3 


50 


0.72 


0.71 


0.75 


0.50 


6 


2.2 


40 


0.80 


0.76 


0.63 


0.41 


7 


3.0 


>240 


0.48 


0.89 


0.69 


0.57 


8 


2.2 


40 


0.89 


1.29 


0.77 


0.87 


9 


2.2 


40 


0.77 


0.77 


0.77 


0.33 


10 


2.8 


150 


0.58 


0.79 


0.75 


0.27 


11 


2.1 


30 


0.91 


0.80 


0.74 


0.60 


12 


1.9 


20 


0.86 


0.80 


0.80 


0.75 


13 


2.1 


30 


0.85 


1.11 


0.70 


0.58 


14 


2.9 


150 


0.55 


0.55 


0.71 


<0 


15 


2.6 


100 


0.72 


0.85 


0.75 


0.38 


16 


2.3 


50 


0.73 


0.78 


0.77 


0.53 



where N is the linear dimension of the image in pixels, and 
the wavelet spectrum as 



(10) 



E{a)= I^K*)! 2 - 

This spectrum, which can be thought of as a smoothed 
version of the power spectrum, reveals the dominant energy 
scales in an image. 

Finally we define the wavelet cross-correlation coeffi- 
cient as 



r w {a) 



£X;W/i(a,x)VF2(a,x) 



(11) 



[£i(a)£ 2 (a)p/2 

where the subscripts refer to two images of the same size. 
This is the wave let analogue of the Fourier cross power spec- 
trum. Following lFrick et al.M200ll) . we assign an uncertainty 
of 



Ar w (a) 



(12) 



Vn~~2 

to this coefficient, where n = (L/a) 2 and L is the linear size 
of the image. 



3.4 Wavelet Cross-correlations 

The wavelet cross-correlation spectrum was calculated 
for each of the following map pairs: (a) radio and 60 )im 
emission; (b) radio and 100 /im emission; (c) radio and 
CO emission; (d) 60 /jm and CO emission; (e) radio and 
Hi emission; (f) 60 fim and Hi emission; (g) 100 /im 
and Hi emission; (h) CO and Hi emission. The median 
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Figure 7. Wavelet cross-correlation spectra for emission from the whole LMC. Top Left. Median filtered 1.4 GHz radio vs. 60 (im 
emission and median filtered 1.4 GHz radio vs. 100 /^m emission. Top Right. Median filtered 1.4 GHz radio vs. integrated CO emission 
and 60 /im vs. CO emission. Bottom Left. Median filtered 1.4 GHz radio vs. integrated Hi emission and 60 fim vs. Hi emission. Bottom 
Right. 100 fim vs. Hi emission and CO vs. Hi emission. For all panels, the wavelet scale is indicated on the bottom axis, and the 
corresponding spatial scale (for an assumed LMC distance of 50.1 kpc) is shown on the top axis. 



filtered 1.4 GHz map was used for the comparisons in- 
volving radio emission. We calculate the cross-correlation 
spectra on wavelet scales log(a) = 1.8 — 3.8 for pair (e), 
log(a) = 1.9 — 3.8 for pairs (a) and (f), and log(a) = 2.3 — 3.8 
for pairs (b), (c), (d), (g) and (h). For an LMC distance of 
50.1 kpc, the range of spatial scales that we probe with all 
eight cross-correlation spectra is 0.05 — 1.5 kpc, extending 
down to 0.02 kpc for the highest resolution datasets. The 



smallest scale that we can study is limited by the resolution 
of the data. The largest scale depends on the size of the 
field, and is determined such that the width of each map is 
equal to or greater than four times the largest wavelet scale. 
The results for all eight pairs are shown in Fig. |7| for the 
whole LMC. Figs. |H| and [5] show the results for pairs (a), (c), 
and (e) for the individual sub-regions, as well as the wavelet 
cross-correlation between the original (unfiltered) radio 
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image with the 60 /tm emission. For the cross-correlation 
spectra of the sub-regions, the largest wavelet scale that we 
study is log(a) = 3.0, corresponding to a spatial scale of 
~0.25 kpc. 



3.4-1 Cross-correlations between radio and FIR emission 

For the whole LMC, it is clear that the correlation between 
the 1.4 GHz radio emission and the 60 or 100 /xm dust 
emission (Fig. [7| top left) is very strong, and significantly 
better at all spatial scales than t he correlations tha t involve 
the molecular or atomic gas. iFrick et alJ l)200l|) regard 
values higher than 0.75 to be indicative of an excellent 
correlation between two images: by this measure, the 
correlation between the radio and 60 //m emission in the 
LMC remains strong down to spatial scales of ~50 pc (210 
arcsec). The correlation between the 1.4 GHz radio and the 
100 /im dust emission is marginally worse across all scales, 
falling below 0.75 at ~60 pc (250 arcsec). Nonetheless, 
it still appears that, averaged over the entire LMC, the 
radio-FIR correlation is maintained on scales less than the 
scale height of the gaseous disc, which has been estimated 
to be between 80-1 00 pc fclmeereen et al.ll20oih and 180 
pc jKim et al]ll999lL 

Since the emission of the LMC as a whole is strongly 
influenced by the 30 Dor region, we have also compared the 
median filtered 1.4 GHz and the 60 fj,m emission for indi- 
vidual sub-regions across the LMC (Fig. |H] filled squares). 
The spatial scale at which the correlation falls below 0.75 
for each sub-region is listed in Table [3] Generally the 
correlation remains excellent, with half of the sub-regions 
showing a strong correlation down to spatial scales of 
<50 pc (log(a) ~ 2.3). However, in several of the other 
sub-regions the correlation breaks down on much larger 
scales (>100 pc, corresponding to log(a) ~ 2.6). We demon- 
strate in Section 13.51 that there is an inverse relationship 
between the size scale of the correlation breakdown and the 
thermal fraction of the radio emission. Thus, the excellent 
correlation between the radio and 60 fim emission that 
we see in sub-regions such as 8 and 12 may be dominated 
by a correlation between the warm dust and the thermal 
radio emission, which does not drive the correlation be- 
tween the total radio and FIR flux density of spiral galaxies. 

Our results for the wavelet cross-correlations are fairly 
robust whether or not median filtering is used to suppress 
compact point sources. When the original (unfiltered) radio 
map is cross-correlated with the 60 /jm map, the correlation 
generally degrades (Fig.|S] open circles), especially for fields 
at the edge of the LMC where the total emission is low and 
fewer of the compact point sources are presumably intrinsic 
(e.g. sub-regions 3, 9, and 14). However, the qualitative 
trends remain the same, since the compact sources usually 
dominate the emission only on the smallest wavelet scales. 
In a few cases (e.g. sub-regions 6 and 13) the median 
filtering appears to worsen the correlation, by removing 
emission from compact sources that have bright 60 /xm 
counterparts and are thus likely to be H II regions. 

On the other hand, an incorrect assumption about the 



resolution of the input images (e.g., the HIRAS-processed 
IRAS images) could cause the correlation to break down 
on larger scales than it should. This is especially a concern 
in sub-region 8, which covers 30 Doradus. Wavelet power 
spectra for this sub-region show a steep decline below a 
scale of ~250 arcsec, whereas the power spectra for other 
sub-regions tend to steepen at smaller characteristic scales 
(~120 arcsec), slightly larger than our estimate of the 
HIRAS resolution. Using an IRAS image processed with 
the HIRES deconvolution method (kindly provided to us 
by J. Surace at IPAC) results in a much stronger wavelet 
cross-correlation that exceeds 0.75 down to scales of 100 
arcsec (Fig. |SJ dashed line in R8) . HIRES is an alternative 
re-processing of the IRAS survey, which, like HIRAS, aims 
to improve on the nominal resolution of the original IRAS 
survey data. HIRES uses the Maximum Correlation Method 
rather than the Pyramid Maximum Entropy algorithm 
used by HIRES, and achieves a resolution of ~1 arcmin 
for the 60/xm band. The discrepancy between the results 
of the wavelet cross-correlation for the HIRAS and HIRES 
data in sub-regio 8 means that we cannot rule out the 
possibility that the radio-FIR correlation remains strong 
on even smaller scales than we deduce here, especially in 
the 30 Dor region, which constitutes a large fraction of the 
galaxy's emission. Future high-resolution data sets from 
Spitzer should be able to overcame these difficulties with 
the IRAS data. 



3.4-2 Cross-correlations between other wavebands 

Cross-correlations involving the CO and Hi images for the 
whole LMC are shown in the last 3 panels of Fig. [7| For 
each image pair, the value of the wavelet cross-correlation 
at wavelet scales log(a) = 2.0, 2.5, 3.0 and 3.5 (correspond- 
ing to spatial scales of 25, 80, 245 and 770 pc) is listed 
in Table |3] The relative weakness of these correlations 
compared to the radio-FIR correlations is immediately 
apparent. Indeed, the strength of the radio-FIR correlation 
is such that correlating one of the gas tracers with either 
the radio or FIR emission yields basically the same result. 
There is also a strong resemblance between the 100 /im-H I 
and CO-H I correlations which is not so easily explained, 
given the lack of a strong correlation between the 100 ^m 
and CO emission. One possibility is that both the 100 ^m 
and CO emissions reflect the typical sizes of star-forming 
regions, but in different evolutionary states, whereas the H I 
emission is dominated by much larger structures. 

Within individual sub-regions, cross-correlations of the 
median filtered 1.4 GHz with the CO and Hi are shown in 
Fig. [5] As expected, for most regions these correlations are 
much weaker than the cross-correlation of 60 /im with 1.4 
GHz radio. Near 30 Doradus (sub- region 8), the integrated 
Hi emission is weakly anti-correlated (—0.4 < r w < 0) with 
the FIR and radio emission across large and small scales. 
There is also a weak anti-correlation between the Hi and 
radio emission at log(a) ~ 2.7 — 2.8 for the field containing 
Nil (sub-region 13); this localised anti-correlation corre- 
sponds to the dense, clumpy ring of molecular and atomic 
gas that surrounds the bright 60 /im and radio source at 
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Figure 8. Wavelet cross-correlation spectra in individual sub-regions for the 1.4 GHz radio vs. HIRAS 60 /im emission and median 
filtered 1.4 GHz radio vs. HIRAS 60 /im emission. For sub-region 8, we also plot the wavelet cross-correlation spectrum for the 1.4 GHz 
radio vs. HIRES 60 jim emission (crosses). The discrepancy between the cross-correlation results for the HIRAS and HIRES data suggests 
that resolution effects may be responsible for the poor correlation between the HIRAS 60 fim and the 1.4 GHz radio emission on small 
scales. 



the centre of the Nil star-forming region. 



3.5 Decomposition of the Radio Continuum 
Emission 

Previous studies of spiral galaxies have shown that the 
local FIR-radio correlation is constituted by a correlation 
between the warm dust and the thermal radio emission, and 
between the cool dust and the non-thermal radio emission 



llHoernes et al.lll99ct) . In order to investigate whether the 
good correlation on small spatial scales that we found for 
some sub-regions is dominated by a correlation between 
the warm dust and thermal radio emission, we used Parkes 
maps of the LMC at 1.4 and 4.8 GHz to estimate the 
thermal fraction of the radio emission across the LMC. 

First, the single-dish 1.4 and 4.8 GHz radio maps 
were smoothed to a common resolution of 15^2. We then 
estimated the thermal fraction of the radio emission in 
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Figure 9. Wavelet cross-correlation spectra in individual sub-regions for the integrated Hi vs. 1.4 GHz emission and integrated CO vs. 
1.4 GHz emission. The dotted lines show the cross-correlation spectra for the median filtered 1.4 GHz radio vs. 60 lira emission from 
Fig. [5] Note that sub-region 1 is poorly covered by our CO map. 



each sub-reg i on by applying equation 1 from the paper by 
Nikl as et alJ lll997l) to the total flux density of each sub- 
region at 1.4 and 4.8 GHz, 

5 4 .8 , /^4.8\" Qth , t /^.SV""*" /io\ 
p = Jth + Jnth ■ (i-O) 

Here, f t h is the thermal fraction of the radio emission at 
1.4 GHz, f nt h = 1 — fth is the non-thermal fraction, and 
5*1.4 and 54.8 are the flux desnsities of each sub-region at 1.4 
and 4.8 GHz respectively. For the thermal spectral index, 
we adopted Qth =0-1; for the non-th ermal spectral index, we 
used a n th=0.7 tHavnes et al.lll99ll) . This simple method of 



separation assumes that the emission is optically thin, and 
that the synchrotron emission can be described by a single 
power-law over the 1.4 to 4.8 GHz frequency range. We note 
that our assumption of a constant non-thermal spectral 
index may introduce a bias in the estimated thermal 
fraction, since the spectral index is expected to steepen 
as cosmic rays propagate from their acceleration sites. 
If these acceleration sites are near star- forming regions, 
then such regions will have a flatter non-thermal spectral 
indices than we assume, and thus the thermal fraction will 
be overestimated. Similarly, the thermal fraction will be 
underestimated in regions that are far from acceleration 
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Table 4. Pearson's correlation coefficient, r p , and wavelet corre- 
lation coefficient, r w , at 4 different spatial scales between image 
pairs for the whole LMC. The wavelet scales shown are log(a) = 
2.0, 2.5, 3.0 and 3.5, corresponding to spatial scales d = 25, 80, 
245 and 770 pc. See also Figure [3 
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sites. However, the uniform distribution of thermal fractions 
observed across the sixteen sub-regions (Fig. I1UI suggests 
that this bias is not severe on the scales that we analyse here. 

For the radio emission across the whole LMC, we 
derive a thermal fraction of 0.45. This increases to 0.57 if 
we restrict our field-of-view to the area of the galaxy that 
belongs to the 16 sub-regions, suggesting that the diffuse 
radio emission around the outer edges of the LMC on 
average has a higher non-thermal fraction than the radio 
emission within the galaxy. The difference in flux density 
between the 1.4 GHz emission from the 16 sub-regions and 
the whole LMC is 62 Jy, so while the radio emission from 
the galaxy's edges is intrinsically weak, it still makes some 
contribution to the LMC's total flux density (471 Jy). The 
results of our decomposition for individual sub-regions are 
listed in Table |!|] The thermal fraction varies from ~0.9 for 
the field containing 30 Doradus, to ~0.1 in the south-west 
of the LMC (sub-region 2). For sub-regions 1 and 14, our 
decomposition produces an unphysical result, i.e. thermal 
fractions less than zero. These anomalous values for f t h 
arise because we fix a n th=0.7 in our decomposition; both 
fields are dominated by a bright compact radio point source 
for which a steeper non-thermal spectral index may be 
more appropriate. The median thermal fraction for the 
remaining sub-regions is 0.45, which is highe r tha n the 
typical value for spiral galaxies f~0.1. ICondonl fl992). but 
consisten t with radio contin u um observations of other dwarf 
galaxies llKlein et alJ Il99lt ISkillman fc Kleml ll98Sl) . The 
sub-regions with the highest 1.4 GHz flux density generally 
have higher thermal fractions. One exception is sub-region 
4, the field south of 30 Doradus. Here the thermal fraction 
is low (0.26), but the 1.4 GHz emission is relatively strong. 
Previous studies of the polarisation of the radio emission in 
the LMC have shown that the ma g netic field is well-or dered 
in this region jKlein et al.lll993ft . iKlein et alJ J! 989) have 
argued that the higher non-thermal fraction in this part 
of the LMC arises through a tidal interaction between the 
LMC and SMC, which drags the magnetic field of the LMC 
towards the SMC and facilitates the diffusion of relativistic 
particles into the intergalactic medium. 

In the bottom panel of Fig. 1101 we plot the spatial scale 
at which the radio-FIR wavelet cross-correlation falls below 
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Figure 10. Top. The WLS slope of the radio vs. 60 /im rela- 
tion, plotted against the thermal fraction at 1.4 GHz. Sub-regions 
where our decomposition method produces thermal fractions less 
than are plotted as = with open squares. Bottom. The 
wavelet scale at which the radio-FIR wavelet cross-correlation 
falls below 0.75, plotted against the thermal fraction of the ra- 
dio emission at 1.4 GHz. The open circle indicates the result for 
sub-region 8 using the HIRES 60/xm data. 



0.75 against the thermal fraction of each sub-region. For 
sub-region 8, where the decorrelation at small scales may be 
partially due to the resolution of the HIRAS FIR data, we 
also plot the result for the wavelet cross-correlation between 
the 1.4 GHz and HIRES FIR data (open circle). There is 
a clear inverse relationship between these two parameters: 
in sub-regions with a high thermal fraction, the radio-FIR 
wavelet cross-correlation holds to much smaller scales. 
We also see a steepening of the WLS fit to the radio-FIR 
correlation with increasing thermal fraction (top panel of 
Fig. 11011 . which is consistent with the increased FIR/radio 
ratio near star-forming regions discussed in Sections 13.11 
and [HI 
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4 DISCUSSION 

The primary result of this study is that the local radio-FIR 
correlation in the LMC is very good: within some regions of 
the LMC, the 60 /im and 1.4 GHz emission remain strongly 
correlated down to spatial scales that correspond to the res- 
olution of our data, i.e. ~20 pc. On almost all scales that we 
investigate, the radio-FIR correlation is markedly stronger 
than the correlation between either the radio or the FIR 
emission with the gas surface density. However, we also find 
that the best correlation between the radio and FIR emis- 
sion occurs in regions where the thermal fraction of the radio 
emission is high. These regions typically host the brightest 
emission regions in the LMC, where a steeper than linear 
correlation is also found. Elsewhere in the LMC, the radio- 
FIR correlation is flatter, and tends to show a poorer cor- 
relation on small spatial scales. In this section, we discuss 
the results of our analysis in relation to previous studies 
of the radio-FIR correlation in the LMC. We also consider 
our results in the context of recent work on the radio-FIR 
correlation in normal spiral galaxies. 

4.1 Comparison with Previous LMC Work 

The striking morphological resemblance between the FIR 
and radio emission has been noted by several previous 
studies, albeit wi th lower re s olutio n data th a n use d in 
this study (e.g. iKlein et al] Il989t IXu et al.l Il992l) . In 
particular, IXu et alT " i^9^r "have studied the correlation 
between the FIR and 6.3 cm radio continuum emission 
in the LMC using pixel-by-pixel techniques. The limiting 
resolution of their data was 4.8 arcmin, allowing them 
to examine the correlation on spatial scales above 70 pc 
(for g!lmc = 50.1 kp c). Our r e sults are generally in good 
agre ement with the IXu et all il992l) results. The values 
that IXu et al derive for the correlation coefficient 

(r p =0.84) between their FIR and 6.3 cm radio map, and 
their mean FIR/1.4 GHz ratio (q mca n = 2.40) are effectively 
identical to our measurement of these quantities (r p =0.86, 
and g m cdian = 2.45). Likewise, we find no radial dependence 
of the FIR/radio ratio. 

Through a statistical analysis of 35 Hn regions, 
IXu et all <ll992l) have found that the total 6.3 cm radio 
emission associated with an H II region in the LMC is more 
centrally concentrated than the total FIR emission, causing 
the FIR/6.3 cm radio ratio to break down on scales below 
70 pc. The authors argued that the central peak in the 
radio emission - which produces a "dip" in the g-value 
- should be attributed to thermal radio emission from 
the ionizing source. At larger radii (~ 70 — 200 pc), they 
find the FIR/radio ratio to be enhanced, corresponding 
to the characteristic s cale of dust h eating around the 
star- forming region. The lXu et all (Il992l) results were based 
on a statistical analysis of H n regions in their data, but our 
g-map clearly confirms their description of a "dip-and-ring" 
morphology for the largest Hn regions in the LMC. In 
Fig. 1111 we present a radial analysis of the radio and 
60 /im emission associated with Nil. Both the radio and 
the infrared emission show a ring-like morphology: the 
ring structure is constituted by both bright, well-known 
star-forming objects and diffuse, filamentary emission at 



r ~ 50-100 pc from a centre position of a=4 h 56 m 43, 
5=-66°29'09"(J2000). The 60 /im emission is clearly more 
extended than the radio emission, however, peaking at 
larger radii and decaying more gradually towards the 
background level. 

Finally, it is worth noting the discrepancy between our 
results for the slope of the local FIR-radio c orrelation (Fig. [3] 
and Fig. and those of IXu et al] dl992^ . An unweighted 
least squares fit to the IXr^tHldTlli99^) data indicated a 
slope of 0.98±0.02 for the whole LMC, consistent with a 
linear correlation. However, by determining both the WLS 
fit and the OLS bisector to our scatterplot of the radio 
and FIR emission in the whole LMC, and by examining 
the emission in different sub-regions across the galaxy, 
we identify a change in the slope of the correlation. This 
change in slope suggests that there are two distinct physical 
mechanisms responsible for the correlation. A steeper 
relation (slope ~ 1.1 — 1.3) is found for the high-intensity 
pixels in the most active star-forming regions, where the 
thermal fraction of the radio emission is relatively high. A 
flatter relation (slope ~ 0.6 — 0.9) applies more generally to 
the diffuse radio and FIR emission across the LMC, and in 
regions with a greater non-thermal contribution to the total 
radio emission (Fig. 2] an d Fig. 1101 . The most probable 
re ason for the dis crepancy between our results and those 
of IXu et al] lll992l) is our ability to sample emission from 
the LMC more finely, due to the higher angular resolution 
of our data. Higher resolution permits a better distinction 
between compact and extended emission than would be 
possible otherwise. 

Our results for the slope of the correlation between the 
radio and FIR e mission are consistent with a cosmic ray di f- 
fusion scenario ]Bicav fc Heloulll99d : iMurphv et all 120051. 
in the sense that while there is an apparent deficit of local 
radio emission across the LMC, the LMC's total radio flux 
density is not unusually low. In other words, a flatter than 
linear slope for the local correlation is not necessarily in- 
consistent with the linear slope of the FIR-radio correlation 
for integrated galaxy flux densities, since the integrated flux 
density measurements used to determine the global corre- 
lation can effectively recoup the synchrotron emission from 
the cosmic ray electrons that have diffused away from their 
production sites. 

4.2 Origin of the global radio-FIR correlation 

4-2.1 Star formation models 

As discussed in ^ conventional models for the radio-FIR 
correlation attribute both sources of luminosity to massive 
star formation. Objections to this explanation are often 
based on the fact that the FIR is a prompt and relatively 
local measure of the ionising output of massive stars, while 
the radio emission results from cosmic ray electron (CRe~) 
energy losses occurring over a much longer timescale and 
a much larger area. On a globally averaged basis and 
in a steady-state situation, however, a strong correlation 
between the radio and FIR emission might still be expected. 
The sensitivity of synchrotron emissivity to magnetic field 
strength is not a serious obstacle given that the globally 
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Figure 11. Radial profile of the 60 um and 1.4 GHz radio emis- 
sion around Nil, centred on a=4 h 56 m 43 s , <5=-66°29'09". The 
surface brightness in an annulus of radius r, around the central 
position for the 60 fim emission (filled squares) and the 1.4 GHz 
radio emission (open circles) is shown. The curves are normalised 
by their raw value at r ~ 80 pc, and the width of each annulus is 
15 pc. 



integrated emission reflects the total energy emitted by an 
electron in its radiating lifetime rather than its instanta- 
neous emitted power. On local scales, such models predict 
that the radio emission should appear as a smeared version 
of the FIR emission, aside from small-scale differences 
due to temporal fluctuations that violate the steady-state 
assumption. Support for this prediction comes from the 
larger radial scale lengths of radio compa red to FIR discs 
l|Bicav fc Heloulll990l : iMurphv et aT1l2005l) . 

Mod els in this category include the "calorimeter" 
model bv lVolkl lll989t) . in which the production of both UV 
photons and the relativistic electrons is proportional to the 
massive star population. The model assumes that in most 
galaxies the UV photons from massive stars are completely 
absorbed and reprocessed by interstellar dust, and that 
relativistic electrons are trapped in the galaxy and lose 
all their energy through synchrotron or inverse Compton 
scattering. The energy densities of the magnetic field and 
the interstellar radiation field are further assumed to be 
proportional, maintaining a constant ratio of synchrotron 
to inverse Compton losses. A sub sequent, more g enera l 
version of the calorimeter model by iLisenfeld et al.l <ll996t) 
allows for finite escape probabilities for relativistic electrons 
and variable optical depth to UV photons. 

With regards to the LMC, three results seem especially 
relevant. First, the existence of a strong FIR-radio correla- 
tion on local scales appears to contradict the calorimeter 
model, although Fig. 1101 implies that this strong correlation 
is mainly driven by the correlation between thermal radio 
and FIR emission, for which a local correlation is naturally 
expected. Second, the high thermal fraction of radio emis- 
sion, coupled with a higher FIR/radio ratio compared to 
spiral galaxies, suggests that a large fraction of relativistic 
electrons can escape the LMC, violating the optically thick 
assumption. In addition, the LMC should be less optically 
thick to UV photons than ma ssive spiral galaxi es because 
of the lower dust-to-gas ratio llFitzpatricklll986l) and high 
porosity implied by the large number of superbubbles and 



supershells. 

The third argument against the calorimeter model in 
the LMC is the relatively flat spectral index of a nt h=0.7, 
whereas in the case of efficient cosmic ray trapping one 
would expect a nt h ~ 1.1 iNiklas fc Beckll997l) . On the other 
hand, the integrated spectral index is a poor diagnostic of 
CRe - escape if either the contribution of SNRs to the total 
radio emis sion or convective transp ort of cosmic rays is not 
negligible iLisenfeld fc Volld [2000 ) . Additional constraints 
could be provided by examining the variation in spectral 
index with distance from star-forming regions. Clearly, a 
more detailed decomposition of the LMC's radio emission 
into thermal, diffuse synchrotron, and SNR emission, using 
mult i- frequency high-resolutions data sets, is required to 
fully address this issue. 



4-2.2 Magnetic field - gas density coupling models 

An alternative "op t ically thin model" has been proposed 
by iHelou fc Bicavl <ll993l) . Once again, the dust-heating 
UV photons and the cosmic-rays are assumed to have a 
common origin in massive star formation, but in this case 
both photons and cosmic rays can escape the galaxy. The 
radio-FIR correlation is enforced by maintaining equal 
escape rates for photons and CRe - 's, i.e. to oc t CS c/tsync, 
where to is the UV optical depth, t csc is the CRe - escape 
time sca le, and f svn c is t he sy nchrotron loss time scale. In 
practice. IHelou fc Bicavl (1993) achieve this by coupling the 
magnetic field B (important for the synchrotron emissivity) 
to the gas density p (important for UV absorption) by 
the relation B oc p ' 5 , assuming that the typical distance 
that a CRe - travels before escape (£ 6 sc) depends only 
weakly on the disc scale height (ho), l CS c oc ho' 5 . Additional 
assumptions include a nearly con stant dust-to-gas r atio, in 
agreement with observations (e.g. JXu fc Heloull996f) . and a 
proportionality between synchrotron and inverse Compton 
losses, as in the calorimeter model. The physical basis for 
the B-p coupling could be amplificati on of MHD turbule nce 
until energy equipartition is reached llGroves et al]l2003ft . 

Other models in this class differ slightly in their 
prescript ion for cosmic-r a y con finement (and hence number 
density 1 ). IHelou fc Bicavl lll993l) assu me a confinement pro - 
cess such that lesch^ ' 5 is constant, iNiklas fc Beckl jl997l) 
assume equipartition between the CR e - density (ucr) and 
magnetic field, ucr oc B 2 oc p, and iHoernes iet_aL| Jl998') 
assum e a nearly constant CRe - density, as do lMurgiaet^rtl 
(2005). In spite of these differences, which may simply 
reflect our poor understanding of how cosmic rays are 
distributed, the basic condition of B oc p 05 is common 
to all of these models, and is motivated by the sensitive 
dependence of synchrotron emission on the magnetic field 
strength, bearing in mind the large range of magnetic 
field strengths obser ved in galaxies obeying the correlation 
llCondon et alJll99ll) . These models have the advantage of 
naturally producing a local radio-FIR correlation, breaking 
down only on scales ~ l csc due to cosmic ray diffusion. 
The fundamental B-p coupling, moreover, could be valid 
down to the scales of typical gas clouds (~100 pc). Most 
significantly, the coupling can even produce a radio-FIR 
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correlation in cases where the FIR is dominated by cool 
dust heated by older stellar populations, although the 
correlation would not necessarily be linear l|Hoernes et alJ 
Il998fi . 

Do our LMC results support the idea of a coupling 
between magnetic field and gas density? Although there 
is a general morphological correspondence between the 
large-scale features of the neutral gas and the radio and FIR 
maps, the Hi and CO emission are poorly correlated with 
the radio emission on scales below ~1 kpc (Figs. [7| and [5J|. 
Even in regions with significant non-thermal radio fractions, 
the radio-FIR correlation cannot be simply derived from 
a correlation between the Hi (which dominates the total 
gas content) and the radio continuum. Note however that 
the synchrotron emission depends on the number density of 
relativistic elec t rons a s well as the magnetic field strength; 
iHelou fc Bicavl l|l993f) constrain ncn by assuming that 
l CBC oc hg' 5 , but provide only qualitative arguments to justify 
this assumption. If the escape length for cosmic rays in the 
LMC is significantly shorter than a nominal value of ~ 2 
kpc - for example because of a high porosity of the ISM - 
then the galaxy cannot effectively extract energy from the 
cosmic rays it generates, even if B attains values typical of 
spiral galaxies. Thus, it is still possible that B-p coupling 
holds in the LMC, but that the CRe~'s are not able to fully 
explore the dense gas distribution. Direct measurements of 
the local magnetic field stre ngth throughout the L MC, e.g. 
through Faraday rotation JCaensler et al]l2005h . will be 
required to rigorously test the B-p coupling model. 

We can, howe ver, address the suggestion by 
iNiklas fe Beckl il997T) that the FIR emission is cou- 
pled to the gas density via the Kennicutt-Schmidt law, 
Psfr oc (pgas) 1 ' 4 , and that this combined with B-p g!ls 
coupling leads to a tight radio-FIR correlation. Our wavelet 
cross-correlation analysis shows that the FIR emission 
is much more poorly correlated with the gas column 
density than with the radio emission, even in sub-regions 
where the thermal radio fraction is relatively small. Thus, 
it appears unlikely that the Kennicutt-Schmidt law can 
serve as a basis for understanding the radio-FIR correlation. 



4-2.3 Connection of FIR-radio correlation to molecular 
gas 

Recently, iMurgia et al.l (120051) have demonstrated that the 
CO-radio correlation within spiral galaxies is as strong 
as the FIR-radio correlation down to scales of ~100 pc. 
They explain both by arguing that the gas density and 
magnetic field strength are linked via equipartition between 
turbulent and magnetic pressure in a hydrostatic disc. 
However, our results indicate a relatively poor correlation 
between CO and 1.4 GHz radio emission within the LMC 
(Fig. [7J. We identify several possible rea sons for the 
differe nce between our results and those of IMurgia et alJ 
(2005). First, the thermal radio fraction f t h is much higher 
for the LMC than in spiral galaxies, so one expects the 
CO-radio correlation on small scales to behave more like 
the CO-60 fim correlation, breaking down near Hn regions 
due to the photodestruction of molecular clouds. Second, 



the fraction of neutral gas mass in molecular form, / mo i, 
is much smaller in the LMC, ~0.05 assuming a Galactic 
CO-to-H2 conversion factor, compared to >0.5 in the 
centra l regi ons of massive galax ies studied bv lMurgia et alJ 
J2005I) (cf. IWong fe Blit3 120021) . Since CO cannot survive 
without self-shielding, this results in a much lower volume 
filling factor for CO emission, and thus no CO counterpart 
to the extended non-thermal radio emission. Finally, the 
combination of strong external pressures exerted on the 
LMC due to its interaction with the Milky Way halo and 
variable internal pressures due to the expansion of bubbles 
and supershells suggests that simple estimates of the 
hydrostatic pressure will not predict the distribution of CO 
clouds. Indeed, observa tions of the LMC 4 supershell by 
lYamaguchi et ail (I2001T) provide evidence for the formation 
of molecular clouds in the swept-up Hi surrounding the 
supershell. 



5 CONCLUSIONS 

We have examined the correlation between the 1.4 GHz 
radio continuum and 60 pm emission in the LMC, using 
pixel-by-pixel techniques and a wavelet cross-correlation 
method. The correlations between the radio continuum, 
FIR emission and the integrated 12 CO and Hi emission 
have also been studied, in order to assess the relevance 
of gas density for the local radio-FIR correlation. The 
high angular resolution datasets that are available for 
the LMC allow us to probe the radio-FIR correlation on 
spatial scales above ~ 20 pc, and to identify variations 
of the FIR/radio ratio for different interstellar environments. 

The main results of our study are: 

1. Pearson's correlation coefficient shows that the 
1.4 GHz and 60/ 100pm are positively correlated across the 
LMC. The correlation coefficient for 1.4 GHz and 60pm 
emission in the whole LMC is ~ 0.86, indicating a strong 
correlation. The correlation improves in regions where 
there are clear signs of star-forming activity, e.g. near 
30 Doradus. The correlation between the radio and FIR 
emission is significantly stronger than correlations between 
these wavebands and cold gas tracers. 

2. Pixel-by-pixel scatterplots indicate that locally there 
are two independent correlations between the radio and 
FIR emission, both of which are non-linear. A steeper 
than linear correlation (slope ~ 1.1 — 1.3) is found for 
the high-intensity pixels in the most active star-forming 
regions, where the thermal fraction of the radio emission 
is relatively high. A flatter correlation (slope ~ 0.6 — 0.9) 
applies more generally to the diffuse radio and FIR emission 
across the LMC, and in regions with a greater non-thermal 
contribution to the total radio emission. 

3. Both the integrated FIR/radio ratio (q = 2.55) and 
the median FIR/radio ratio across the LMC (^med ian = 2.45) 
are hi gher than the mean FIR/radio ratio of the lYun et alJ 
(2001) galaxy sample (q mC an = 2.34), suggesting that the 
radio emission in the LMC is slightly underluminous, or 
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that there is a slight excess of FIR emission. In contrast 
to nearby spiral galaxies, the brightest structures in 
the FIR/radio ratio map of the LMC do not correspond 
to the brightest structures in the input FIR and radio maps. 

4. A simple decomposition of the 1.4 GHz radio 
emission indicates that the thermal fraction of the LMC i s 
higher than for normal disc galaxies (~ 0.1, ICondonlll99l . 
For an assumed non-thermal spectral index of a nt h=0-7, 
the thermal fraction for the whole the LMC is 0.45, varying 
between and 0.9 for individual sub-regions. The high 
FIR/radio ratio and high thermal fraction of the LMC's 
radio emission suggests that a large proportion of CRe~s 
may be escaping from the LMC with low synchrotron losses. 

5. Wavelet cross-correlation spectra show that the 
correlation between the 1.4 GHz radio continuum and 
60 fim emission is very good across the LMC, even on scales 
corresponding to a few tens of parsecs. For the whole LMC, 
the correlation between the radio and 60 /im emission on 
all scales is significantly better than the correlation of the 
Hi or CO with either the radio or FIR emission. We argue, 
however, that the excellent correlation between the radio 
and FIR emission for the LMC as a whole is driven by the 
30 Doradus region. For individual sub-regions within the 
LMC, there is an inverse relationship between the scale 
on which the correlation breaks down and the thermal 
fraction of the radio emission: regions that show a strong 
correlation to small spatial scales are also the regions where 
the thermal fraction of the radio emission is high. In regions 
where the thermal fraction is less than 0.5, the wavelet 
cross-correlation between radio and 60 fim emission breaks 
down across a range of scales between 50 and 250 pc. 

While models for the radio-FIR correlation typically 
address the emission and coupling mechanisms in normal 
spiral galaxies, the LMC remains an interesting testbed 
because of the wide range of angular scales that can be 
probed. The presence of a strong local radio-FIR correlation 
that is dominated by high-brightness thermal emission, 
combined with evidence for a diffuse non-thermal radio 
component that contributes to the LMC's total radio 
flux density, underscores that different phenomena can be 
responsible for brightness as opposed to integrated flux 
density. This distinction should be kept in mind in future 
studies of more massive galaxies. 
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